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Y. “ Mechanical Integration of the Linear Differential Equations 
of the Second Order with Variable Coefficients.” By Prof. 
Sir William Thomson, LL.D., F.R.S. Received January 
28 , 1876 . 

Every linear differential equation of the second order may, as is 
known, be reduced to the form 

±(L c h\ =u , .. . (i) 

dx VP dx)’ K 1 

where P is any given function of x. 

On account of the great importance of this equation in mathematical 
physics (vibrations of a non-uniform stretched cord, of a hanging chain, 
of water in a canal of non-uniform breadth and depth, of air in a pipe of 
non-uniform sectional area, conduction of heat along a bar of non-uniform 
section or non-uniform conductivity, Laplace’s differential equation of 
the tides, &c. &c.), I have long endeavoured to obtain a means of facili¬ 
tating its practical solution. 

Methods of calculation such as those used by Laplace himself are 
exceedingly valuable, but are very laborious, too laborious unless a serious 
object is to be attained by calculating out results with minute accuracy. 
A ready means of obtaining approximate results which shall show the 
general character of the solutions, such as those so well worked out by 
Sturm*, has always seemed to me a desideratum. Therefore I have made 
many attempts to plan a mechanical integrator which should give solu¬ 
tions by successive approximations. This is clearly done now, when we 
have the instrument for calculating \p(x)dx, founded on my brother’s 
disk-, globe-, and cylinder-integrator, and described .in a previous com¬ 
munication to the Royal Society; for it is easily proved t that if 

H 1 dx)dx, ] 

r l .(2) 

=J* P(C-£ U 2 dx) dx, J 

&c., 

where u x is any function of x, to begin with, as for example u x ~x• then 
u 2 , u 8 , &c. are successive approximations converging to that one of the 
solutions of (1) which vanishes when #=0. 

Now let my brother’s integrator be applied to find C—J 0 * u. L dx , and let its 
result feed, as it were, continuously a second machine, which shall find 
the integral of the product of its result into P dx. The second machine 

* “ Memoire sur les equations differentielles lineaires dip second ordre,” Liouville’s 
Journal, vol. i. 1836. 

t Cambridge Senate-House Examination, Thursday afternoon, January 22nd, 1874. 
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will give out continuously the value of u 2 . Use again the same process 
with u 2 instead of u v and then w 3 , and so on. 

After thus altering, as it were, u x into u 2 by passing it through the 
machine, then u 2 into u 3 by a second passage through the machine, and so 
on, the thing will, as it were, become refined into a solution which will 
be more and more nearly rigorously correct the oftener we pass it through 
the machine. If u i+1 does not sensibly differ from u i} then each is sen¬ 
sibly a solution. 

So far I had gone and was satisfied, feeling I had done what I wished 
to do for many years. But then came a pleasing surprise. Compel 
agreement between the function fed into the double machine and that 
given out by it. This is to be done by establishing a connexion which 
shall cause the motion of the centre of the globe of the first integrator of 
the double machine to be the same as that of the surface of the second 
integrator’s cylinder. The motion of each will thus be necessarily a 
solution of (1). Thus I was led to a conclusion which was quite unex¬ 
pected ; and it seems to me very remarkable that the general differential 
equation of the second order with variable coefficients may be rigorously, 
continuously, and in a single process solved by a machine. 

Take up the whole matter ah initio : here it is. Take two of my 
brother’s disk-, globe-, and cylinder-integrators, and connect the fork 
which guides the motion of the globe of each of the integrators, by 
proper mechanical means, with the circumference of the other inte¬ 
grator’s cylinder. Then move one integrator’s disk through an angle =x, 
and simultaneously move the other integrator’s disk through an angle 
always =j* 0 *P dx, a given function of x. The circumference of the 
second integrator’s cylinder and the centre of the first integrator’s globe 
move each of them through a space which satisfies the differential 
equation (1). 

To prove this, let at any time g v g 2 be the displacements of the centres 
of the two globes from the axial lines of the disks ; and let dx, P dx be 
infinitesimal angles turned through by the two disks. The infinitesimal 
motions produced in the circumferences of two cylinders will be 

g x dx and </ 2 P dx. 

But the connexions pull the second and first globes through spaces 
respectively equal to those moved through by the circumferences of the 
first and second cylinders. Hence 

g x dx—dg 2 t and gfddx~dg 1 ; 

and eliminating g 2 , 

±(±<* 9 ? W/ 

dx\Sdco) Jv 

which shows that g x put for u satisfies the differential equation (1). 
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The machine gives the complete integral of the equation with its two 
arbitrary constants. Tor, for any particular value of a?, give arbitrary 
values Gq, Gt 2 . [That is to say mechanically; disconnect the forks from the 
cylinders, shift the forks till the globes’ centres are at distances Gq, Gr a 
from the axial lines, then connect, and move the machine.] 

We have for this value of oc, 

9i= a j) andg=G 2 P; 


that is, we secure arbitrary values for g 1 and by the arbitrariness of 

ctoo 

the two initial positions Gq, Gq of the globes. 


YI. “ Mechanical Integration of the general Linear Differential 
Equation of any Order with Variable Coefficients/'’ By Prof. 
Sir William Thomson^ LL.D., T.R.S. Received January 2S 3 
1876. 

Take any number i of my brother’s disk-, globe-, and cylinder-inte¬ 
grators, and make an integrating chain of them thus :—Connect the 
cylinder of the first so as to give a motion equal to its own* to the fork 
of the second. Similarly connect the cylinder of the second with the 
fork of the third, and so on. Let g v g 2 , g 3 , up to g i9 be the positions t of 
the globes at any time. Let infinitesimal motions P flat, P flee, P 3 dx ,.... 
be given simultaneously to all the disks (doc denoting an infinitesimal 
motion of some part of the mechanism whose displacement it is convenient 
to take as independent variable). The motions (dtc 19 dic 2 ,... cfo*) of the 
cylinders thus produced are 

dKy^g^jdati ^ 2 =^PA • * • dKi—g^flx. , . . (1) 

But, by the connexions between the cylinders and forks which move the 
globes, dK x ~dg 29 dic 2 = dg 3 ,,,. d^^dg^, and therefore 

dg^cji^dx, dg 3 =g 2 ¥ 2 dx,...clg i =cj i _ l ¥ i _ 1 dx ) 

[ ■ • ( 2 ) 

dK 1 —g 1 P^, dK 2 -=g 2 ~P 2 doo, ... dK i ~g i l y i dx. J 
n _1 d I d 1 d 1 dXi 

(Jl ^TgTx¥jh!'‘'¥~ x dx¥ i d^ .^ ’ 

Suppose, now, for the moment that we couple the last cylinder with the 

* For brevity, the motion of the circumference of the cylinder is called the cylinder’s 
motion. 

t For brevity, the term “ position ” of any one of the globes is used to denote its 
distance, positive or negative, from the axial line of the rotating disk on which it 


and 

Hence 




